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Abstract 

Wightman function and vacuum expectation value of the field square are evaluated for 
a massive scalar field with general curvature coupling parameter subject to Robin boundary 
conditions on two codimension one parallel branes located on (Z?+l)-dimcnsional background 
spacetime AdSut+i x E with a warped internal space E. The general case of different Robin 
coefficients on separate branes is considered. The application of the generalized Abel-Plana 
formula for the series over zeros of combinations of cylinder functions allows us to extract 
manifestly the part due to the bulk without boundaries. Unlike to the purely AdS bulk, 
the vacuum expectation value of the field square induced by a single brane, in addition to 
the distance from the brane, depends also on the position of the brane in the bulk. The 
brane induced part in this expectation value vanishes when the brane position tends to the 
AdS horizon or AdS boundary. The asymptotic behavior of the vacuum densities near the 
branes and at large distances is investigated. The contribution of Kaluza-Klein modes along 
£ is discussed in various limiting cases. In the limit when the curvature radius for the AdS 
spacetime tends to infinity, we derive the results for two parallel Robin plates on background 
spacetime R( Dl,1 > x E. For strong gravitational fields corresponding to large values of the 
AdS energy scale, the both single brane and interference parts of the expectation values 
integrated over the internal space are exponentially suppressed. As an example the case 
E = S 1 is considered, corresponding to the AdSo+i bulk with one compactified dimension. 
An application to the higher dimensional generalization of the Randall-Sundrum brane model 
with arbitrary mass terms on the branes is discussed. 



PACS numbers: 04.62. +v, ll.10.Kk, 04.50.+h 



1 Introduction 

Anti-de Sitter (AdS) spacetime is remarkable from different points of view. The early interest 
to this spacetime was motivated by the questions of principal nature related to the quantization 
of fields propagating on curved backgrounds. The investigation of the dynamics of fields on 
AdS is of interest not only because AdS space is a space of high symmetry, and hence exact 
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solutions for the free field theory can be written down, but also because it is a space of con- 
stant nonzero curvature, and thus field theory in AdS background is not a trivial rewriting of 
Minkowski spacetime field theory. The presence of the both regular and irregular modes and 
the possibility of interesting causal structure lead to many new phenomena. The importance 
of this theoretical work increased when it has been discovered that AdS spacetime generically 
arises as ground state in extended supergravity and in string theories. Further interest in this 
subject was generated by the appearance of two models where AdS geometry plays a special 
role. The first model, the so called AdS/CFT correspondence (for a review see P), repre- 
sents a realization of the holographic principle and relates string theories or supergravity in the 
AdS bulk with a conformal field theory living on its boundary. It has many interesting con- 
sequences and provides a powerful tool to investigate gauge field theories. The second model 
is a realization of a braneworld scenario with large extra dimensions and provides a solution 
to the hierarchy problem between the gravitational and electroweak mass scales (for reviews in 
braneworld gravity and cosmology see Ref. |2J). In this model the main idea to resolve the 
large hierarchy is that the small coupling of four dimensional gravity is generated by the large 
physical volume of extra dimensions. Braneworlds naturally appear in string/M-theory context 
and provide a novel setting for discussing phenomenological and cosmological issues related to 
extra dimensions. The model introduced by Randall and Sundrum 3 is particularly attractive. 
The corresponding background solution consists of two parallel flat branes, one with positive 
tension and another with negative tension embedded in a five dimensional AdS bulk. The fifth 
coordinate is compactified on S 1 /Z2, and the branes are on the two fixed points. In the original 
version of the model it is assumed that all matter fields are confined on the branes and only the 
gravity propagates freely in the five dimensional bulk. More recently, alternatives to confining 
particles on the brane have been investigated and scenarios with additional bulk fields have been 
considered. Apart from the hierarchy problem it has been also tried to solve the cosmological 
constant problem within the braneworld scenario. The braneworld theories may give some al- 
ternative discussion of the cosmological constant. The basic new ingredient is that the vacuum 
energy generated by quantum fluctuations of fields living on the brane may not curve the brane 
itself but instead the space transverse to it. 

Randall- Sundrum scenario is just the simplest possibility within a more general class of higher 
dimensional warped geometries. Such a generalization is of importance from the viewpoint of 
a underlying fundamental theory in higher dimensions such as ten dimensional superstring the- 
ory. From a phenomenological point of view, higher dimensional theories with curved internal 
manifolds offer a richer geometrical and topological structure. The consideration of more gen- 
eral spacetimes may provide interesting extensions of the Randall- Sundrum mechanism for the 
geometric origin of the hierarchy. Spacetimes with more than one extra dimension can allow for 
solutions with more appealing features, particularly in spacetimes where the curvature of the 
internal space is non-zero. More extra dimensions also relax the fine-tunings of the fundamental 
parameters. These models can provide a framework in the context of which the stabilization 
of the radion field naturally takes place. Several variants of the Randall-Sundrum scenario in- 
volving cosmic strings and other global defects of various codimensions has been investigated in 
higher dimensions (see, for instance, [I]- JE] and references therein). In particular, much work 
has been devoted to warped geometries in six dimensions (see references in |17j). 

In braneworld models the investigation of quantum effects is of considerable phenomeno- 
logical interest, both in particle physics and in cosmology. The braneworld corresponds to a 
manifold with boundaries and all fields which propagate in the bulk will give Casimir-type con- 
tributions to the vacuum energy (for reviews of the Casimir effect see Ref. [18]). and as a result 
to the vacuum forces acting on the branes. Casimir forces provide a natural alternative to the 
Goldberger-Wise mechanism for stabilizing the radion field in the Randall-Sundrum model, as 
required for a complete solution of the hierarchy problem. In addition, the Casimir energy gives 
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a contribution to both the brane and bulk cosmological constants and, hence, has to be taken 
into account in the self-consistent formulation of the braneworld dynamics. Motivated by these, 
the role of quantum effects in braneworld scenarios has received a great deal of attention |19j - 
|49j . However, in the most part of these papers the authors consider the global quantities such 
as the total Casimir energy, effective action or conformally invariant fields. The investigation 
of local physical characteristics in the Casimir effect, such as expectation value of the energy- 
momentum tensor and the field square, is of considerable interest. Local quantities contain 
more information on the vacuum fluctuations than the global ones. In addition to describing 
the physical structure of the quantum field at a given point, the energy-momentum tensor acts 
as the source in the Einstein equations and therefore plays an important role in modelling a 
self-consistent dynamics involving the gravitational field. Quantum fluctuations play also an 
important role in inflationary cosmology as they are related to the power spectrum. As it has 
been shown in Ref. jlS]) the quantum fluctuations of a bulk scalar field coupled to a brane 
located scalar field with a bi-quadratic interaction generate an effective potential for the field on 
the brane with a true vacuum at the nonzero values of the field. In particular, these calculations 
are relevant to the bulk inflaton model (SHI) where the inflation on the brane is driven by the 
bulk scalar field. In the case of two parallel branes on AdS background, the vacuum expectation 
value of the bulk energy-momentum tensor for a scalar field with an arbitrary curvature coupling 
is investigated in Refs. jSHEHl- In particular, in Ref. the application of the generalized 
Abel-Plana formula |51[ 152} to the corresponding mode sums allowed us to extract manifestly 
the parts due to the AdS spacetime without boundaries and to present the boundary induced 
parts in terms of exponentially convergent integrals for the points away the boundaries. The 
interaction forces between the branes are investigated as well. Depending on the coefficients in 
the boundary conditions, these forces can be either attractive or repulsive. The local vacuum 
effects for a bulk scalar field in brane models with dS branes are discussed in Refs. |4fi| 14=7] . 
On background of manifolds with boundaries, the physical quantities, in general, will receive 
both volume and surface contributions. For scalar fields with general curvature coupling, in 
Ref. ]53] it has been shown that in the discussion of the relation between the mode sum energy, 
evaluated as the sum of the zero-point energies for each normal mode of frequency, and the 
volume integral of the renormalized energy density for the Robin parallel plates geometry it is 
necessary to include in the energy a surface term located on the boundary. An expression for the 
surface energy-momentum tensor for a scalar field with general curvature coupling parameter 
in the general case of bulk and boundary geometries is derived in Ref. [3]. The corresponding 
vacuum expectation values are investigated in Ref. 05] for the model with two flat parallel 
branes on AdS bulk and in 0f5] for a dS brane in the flat bulk. In particular, for the first case 
it has been shown that in the Randall- Sundrum model the cosmological constant induced on 
the visible brane by the presence of the hidden brane is of the right order of magnitude with 
the value suggested by the cosmological observations without an additional fine tuning of the 
parameters. 

The purpose of the present paper is to study the Wightman function and the vacuum expec- 
tation value of the field square for a scalar field with an arbitrary curvature coupling parameter 
obeying Robin boundary conditions on two codimension one parallel branes embedded in the 
background spacetime AdSpj+i X £ with a warped internal space E. The quantum effective 
potential and the problem of moduli stabilization in the orbifolded version of this model with 
zero mass parameters on the branes are discussed recently in Ref. 0^. In particular, it has been 
shown that one loop-effects induced by bulk scalar fields generate a suitable effective potential 
which can stabilize the hierarchy without fine tuning. The corresponding results are extended to 
the type of models with unwarped internal space E in Ref. 02] ■ We organize the present paper 
as follows. In the next section we evaluate the Wightman function in the region between the 
branes. By using the generalized Abel-Plana formula, we present this function in the form of 
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a sum of the Wightman function for the bulk without boundaries and boundary induced parts. 
The vacuum expectation value of the field square for a general case of the internal space £ is 
discussed in section El for the case of a single brane geometry and in section 0] in the region be- 
tween two branes. Various limiting cases are discussed when the general formulae are simplified. 
A simple example with the internal space S 1 is then considered in section [21 The last section 
contains a summary of the work. The vacuum expectation values of the energy-momentum 
tensor and the interaction forces between the branes will be discussed in the forthcoming paper. 

2 Wightman function 

Consider a scalar field <p(x) with curvature coupling parameter £ satisfying the equation of 
motion 

{g MN V M Vn + m 2 + (R) <p(x) = 0, (2.1) 

M, N = 0, 1, . . . , D, with R being the scalar curvature for a (D + l)-dimensional background 
spacetime, V m is the covariant derivative operator associated with the metric tensor gMN (we 
adopt the conventions of Ref. |55j for the metric signature and the curvature tensor). For mini- 
mally and confer mally coupled scalars one has £ = and £ = Cd = (D — 1) /4D correspondingly. 
We will assume that the background spacetime has a topology AdSo 1 +i x S, where £ is a 
/^-dimensional manifold. 

First let us consider a general class of spacetimes described by the line element (see also the 
discussion in Ref. |4*T] ) 

ds 2 = g MN dx M dx N = e-^rj^dx^dx" - e- 2p(v) j ij dX i dX j - dy 2 , (2.2) 

with rj^ u = diag(l, —1, . . . , —1) being the metric for the Di-dimensional Minkowski spacetime 
r( d i~ i' 1 ) and the coordinates X 1 cover the manifold E, D = D% + Z?2- Here and below /i, v = 
0,1, . . . , D\ — 1 and i,j = 1, . . . , D2. The scalar curvature for the metric tensor from (j2.2|) is 
given by the expression 

R = -F(y)-e 2 PR (l) , (2.3) 
where we have introduced the notation 

F(y) = -2r"(y) + r' 2 {y) + D^iy) + D 2 p' 2 (y), (2.4) 

with r(y) = D\a(y) + D2p(y) and R(~\ is the scalar curvature for the metric tensor 7^. The 
background metric for many braneworld models falls into the general class corresponding to Eq. 
(|2.2|) . In particular, the metric with D2 = and o{y) = ko\y\ is the Randall-Sundrum case. 
The metric solutions with D2 = 1, a(y) = kz>\y\, p(y) = const and D% = 1, cr(y) = p(y) = kz>\y\ 
correspond to global LQ and local 5 a strings respectively. In the case of more extra dimensions 
generalizations of these models with S = S Di were considered in [71 1111 113j . The general 
case of a symmetric homogeneous space S with topologically non-trivial Yang-Mills fields is 
investigated in Ej. In particular, when the bulk cosmological constant is negative, a simple 
solution exists with u{y) = p(y) = kpy. For the case of S with nonzero Ricci scalar one can 
have a regular solution in the bulk with p(y) = const and Randall-Sundrum type warp in the 
Minkowski direction. The phenomenology of this type of scenario has been considered in |16j . 

Our main interest in this paper will be the Wightman function and the vacuum expectation 
value (VEV) of the field square brought about by the presence of two parallel infinite branes of 
codimension one, located at y = a and y = b, a < b. We will assume that on this boundaries 
the scalar field obeys the mixed boundary conditions 

(Ay + Bydy^j (p(x) = 0, y = o, 6, (2.5) 
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with constant coefficients A y , B y . In the orbifolded version of the model which corresponds 
to a higher dimensional Randall-Sundrum braneworld these coefficients are expressed in terms 
of the surface mass parameters and the curvature coupling of the scalar field (see below). The 
imposition of boundary conditions on the quantum field modifies the spectrum for the zero-point 
fluctuations and leads to the modification of the VEVs for physical quantities compared with the 
case without boundaries. The resulting quantum effects can either stabilize or destabilize the 
branewolds and have to be taken into account in the self-consistent formulation of the braneworld 
dynamics. As a first stage in the investigations of local quantum effects, in this section we will 
consider the positive frequency Wightman function defined as the expectation value 

G + (x,x>) = {0\ip(xMx')\0), (2.6) 

where |0) is the amplitude for the vacuum state. The VEVs of the field square and the energy- 
momentum tensor can be obtained from the Wightman function in the coincidence limit with 
an additional renormalization procedure. Instead of the Wightman function we could take any 
other two-point function. The reason for our choice of the Wightman function is related to 
the fact that this function also determines the response of particle detectors in a given state of 
motion. Let (p a (x) be a complete set of positive frequency solutions to the field equation (|2.1j) 
satisfying boundary conditions (j2.5|) and a denotes a set of quantum numbers specifying the 
solution. Expanding the field operator over this set of eigenfunctions and using the commutation 
relations, the Wightman function is presented as the mode sum: 

G + (x,x') = Y,Mx)ip* Q (x'). (2.7) 

a 

The symmetry of the bulk and boundary geometries under consideration allows to present 
the corresponding eigenfunctions in the decomposed form 

va{* M ) =M^)M x )fn(y), (2-8) 

with the standard Minkowskian modes in 

k (^) = - =, k? = (w, k), 



uj =yk 2 + m| n , k = |k|. 

The separation constants mp n are determined by the geometry of the internal space S and by 
the boundary conditions imposed on the branes and will be given below. The modes ^(X) are 
eigenfunctions for the operator A( 7 ) + C,R(^)'- 

[A (7) + 0R (7) ] MX) = -A|#(X), (2.10) 
with eigenvalues A| and normalization condition 



J d D *X^ p { X )^U X ) = 5 W' 



(2.11) 



A( 7 ) is the Laplace-Beltrami operator for the metric 7^. In the consideration below we will 
assume that A/3 > 0. Substituting eigenfunctions Q2.8[) into the field equation ()2.1|) and using 
equation (|2.1Uj) . for the function f n (y) one obtains the following equation 

l-e T d y e- T d y + m 2 - (F(y) + A^e 2 "] f n = m% n e 2 ° f n . (2.12) 
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Equation is valid for any warp factors cr(y) and p{y). In this paper we consider the case 

of two equal warp factors, with 

p(y) = <?{y) = k D y- (2.13) 

In the case when £ is a torus, this corresponds to a toroidal compactification of AdSp+i space- 
time. For warp factors ()2.13|> equation (|2.12|) takes the form 

[-eP a aye- Do d y + m 2 - QD(D + 1)4] /„ = (m% n - X 2 p ) e 2fJ f n . (2.14) 

The operator in the left hand side does not depend on the internal quantum numbers (3. Thus 
in this case the eigenfunctions f n and the combination rr? n = m| n — X\ do not depend on 0. As 
a result, the dependence on (3 and n of the masses is factorized (see also [IT]): 

m% n = ml + X\. (2.15) 

For the region between the branes, a < y < b, the solution to equation ()2.14|) satisfying the 
boundary condition on the brane y = a is given by expression 

f n (y) = C n e DcT/2 g L ,(m n z a ,m n z), (2.16) 

where 

g v (u,v) = J v {v)Yy*{u) - Ji a \u)Y v {v), (2.17) 
Jv(x), Y v {x) are the Bessel and Neumann functions, 

pkDV pkD3 

j = a,b, (2.18) 



kD 3 kf) 
and 



\D/2) 2 - D{D + l)( + m 2 /k 2 D . (2.19) 
In formula l)2.17j) for a given function F(x) we use the notation 

F {j) (x) = AjF(x) + BjxF'(x), j = a, b, (2.20) 

with the coefficients 

Aj = Aj + Bjk D D/2, Bj = Bjk D . (2.21) 

In the discussion below we will assume values of the curvature coupling parameter for which 
v is real. For imaginary v the ground state becomes unstable |5fi| . For a conformally coupled 
massless scalar one has v = 1/2 and the cylinder functions in Eq. (|2,17|) are expressed via the 
elementary functions. For a minimally coupled massless scalar v = D/2 and the same is the 
case in odd spatial dimensions. Note that in terms of the coordinate z introduced by relation 
(|2 . 18|) , the metric tensor from (|2,2j) with warp factors (|2.13|) is conformally related to the metric 
of the direct product space R( Dl '^ x S by the conformal factor (k£>z)~ 2 . 

From the boundary condition on the brane y = b we receive that the eigenvalues m n have to 
be solutions to the equation 

C«\z b /z a , m n z a ) = jl a \m n z a )Y^(m n z h ) - Y^ \m n z a )J<f> \m n z b ) = 0. (2.22) 

This equation determines the tower of radial Kaluza-Klein (KK) masses. We denote by z = ^y u ,n, 
n = 1,2,..., the zeros of the function C^ b (r], z) in the right half-plane of the complex variable 
z, arranged in the ascending order, j Vin < 7„ )n+ i. The eigenvalues for m n are related to these 
zeros as 

m n = k D ^ n e~ kDa = j^n/za. (2.23) 
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From the orthonormality condition for the functions f n (y) for the coefficient C n in Eq. (|2.16jl 
one finds 

C 2 n = p^T?( nn ), (2.24) 
2k D z£ z a 

where we have introduced the notation 

T ^,u) = u (45^ l A l + B l^u* - u 2 )] -Al + 5 a V " - 2 )| • (2-25) 

I Jv \J]U) J 

Note that, as we consider the quantization in the region between the branes, z a < z < Zb, the 
modes defined by (|2.1(ij) are normalizable for all real values of v from Eq. (|2.19fl . 

Substituting the eigenfunctions (|2.8() into the mode sum (|2.7|) . for the expectation value of 
the field product one finds 

a /3 

/do 
dke ,k(x-x')^ V(>n)r a fe(r?)>n)) (226) 
n=l 

where x = (x 1 , x 2 , . . . , x Dl_1 ) represents the spatial coordinates in R^ Dl ~ l,1 \ and 



hf3v{u) = ug u (u,uz/z a )g u (u,uz' /z a )—= . (2.27) 

lu 2 /zl + k 2 + \ 2 p 

As the expressions for the zeros j u ^ n are not explicitly known, the form (|2.26|) of the Wightman 
function is inconvenient. In addition, the terms in the sum over n are highly oscillatory functions 
for large values n. It is possible to overcome both these difficulties by applying to the sum over 
n the summation formula derived in Refs. |51| I52j by making use of the generalized Abel-Plana 
formula. For a function h(u) analytic in the right half-plane Heu > this formula has the form 

7r 2 ^ f°° h(x)dx 

yL^^^v) = / (a)2 — 

Z n=l Jo Jv {X) + Y) (x) 



dxQ au (x,r]x) h(xei) + h(xe 2) ; (2.28) 



where 



1UC.)- . M f W/ ^ W , a) ■ (2.29) 
K v a \u)li h \v)-K v h \v)I v a \u) 

I u {u) and K v {u) are the Bessel modified functions. Formula (|2.28|) is valid for functions fo(u) 
satisfying the conditions 

\h(u)\ < ei(x)e Cl|?/l |«|^oo, u = x + iy, (2.30) 

and fo(«) = o(n^ 1 ), -u — > 0, where c\ < 2(r] — 1), x 2SBa0 ~ 1 ei(x) — > for x — > +00. Using 
the asymptotic formulae for the Bessel functions for large arguments when 1/ is fixed (see, e.g., 
we can see that for the function hp v {u) from Eq. (|2.27j) the condition (j2.30j) is satisfied if 
z + z' + \ t — 1'\ < 2zb- In particular, this is the case in the coincidence limit t = t' for the region 
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under consideration, z a < z,z' < z^. Applying to the sum over n in Eq. (|2.26|) formula (|2.28j) . 
one obtains 



(0|^(x)^(a/)|0> 



hp v (u)du 2 

7T 



1 



-jo 4 a)2 (u) + yr» 



(0)2, 



ai UZb) 



u 2 -k 2 - X 2 



x Gi a \uz a , uz)G^(uz a ,uz') cosh y u 2 -k 2 - X 2 (t-t')^ j. (2.31) 

where we have introduced notations 

G#'>(u, «) = I„(v)K®{u) - lP(u)K u (v), j = a, b, (2.32) 

(the function with j = b will be used below). Note that we have assumed values of the coefficients 
Aj and Bj for which all zeros for Eq. (|2,22|) are real and have omitted the residue terms in the 
original formula in Refs. |51l I52j . In the following we will consider this case only. 

By the way similar to that used in Ref. [221 > ^ can be seen that under the condition 
z + z' — \t — t I > 2z a , the first term in the figure braces in Eq. (|2.31j) is presented in the form 



hp v {u)du 



Z aJ0 j^\ U )+Y v 



r(a)2 i 



U 



f 

Jo 



duu- 



x 2 + k 2 + A| 



J v {uz)J u (uz) 



lj a \uz a ) K v {uz)K v {uz') 
* J^TxJ Ki a) (uz a ) lu 2 -k 2 -X 2 



du u 



2 - k 2 - X 2 Jt - t') 



x cosh 

Substituting this into formula (|2.31|) . for the Wightman function one finds 

<o|v?(*M*')|o) = (ip(xMx'))^ + M*M*')) (a) 

j3 J 
x / duuG^(uz a ,uz)Gl a \uz a ,uz') 



(2.33) 



tt au (uz a ,uz b ) 



cosh 



u 2 - k 2 - X 2 



u 2 - k 2 - X 2 (t - t') 



(2.34) 



Here the term 



(<p(x)(p(x )) 



M\(o) 



k%-\zz') D / 2 



\D 

2^1^!- 



— J2M X WX') J dke* 



k(x-x') 



duu- 



+ k 2 + X 2 



J v {uz)J v {uz) 



(2.35) 
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does not depend on the boundary conditions and is the Wightman function for the AdS£> 1+ i x £ 
spacetime without branes. The second term on the right of Eq. (|2.34j) . 

p 

duu I^ a \uz a ) K u {uz)K v {uz') 



y/*+\% Ki a \uz a ) ^v? — k 2 — \p 



x cosh 



u 2 



k 2 -X 2 (t-t') , (2.36) 



does not depend on the parameters of the brane at z = z b and is induced in the region z > z a by 
a single brane at z = z a when the boundary z = z b is absent. Thus the last term on the right of 
formula ()2.34|) is the part in the Wightman function induced by the presence of the second brane. 
Hence, the application of the summation formula based on the generalized Abel-Plana formula 
allowed us (i) to escape the necessity to know the explicit expressions for the zeros j u ,n, (h) to 
extract from the two-point function the boundary- free and single brane parts, (iii) to present the 
remained part in terms of integrals with the exponential convergence in the coincidence limit. 
By the same way described above for the Wightman function, any other two-point function can 
be evaluated. Note that expression (|2.35f) for the boundary-free Wightman function can also be 
written in the form 



M*Mz')> (0) = k^\zz') D /^MxW^x') 



p 

oo 



x / duuG+ (Dl _ lil) (x >1 ,x"'; ^u 2 + \ 2 )J u {mz)J v {mz), (2.37) 

where G^ Di _ 11 - ) (x fl ,x'^;m) is the Wightman function for a scalar field with mass m in the 
-Di-dimensional Minkowski spacetime. 
By using the identity 

g £ {uZh) I v {uz)I v {uz>) - ^} UZa) K u (uz)K v (uz') = 

n^Qj u (uz a ,uzf ) )G^J' > (uzj,uz)G^\uzj,uz'), (2.38) 

j=a,b 

where n^' = 1 for the region y > j and rv^' = — 1 for the region y < j (hence, in (|2.38f) one has 
n i a ) = i 5 n ( b ) — —i as we consider the region a < y < b), and 

JW«,«) = -n — I 4r^^P n ■ ( 2 - 39 ) 

the Wightman function in the region z a < z < z b can also be presented in the form 

<0K*Ma/)|0) = (<p(x)<p(x , ))V) + (<p(x)<p(x , )f'> 
k°~\zz') D l 2 



p 

Q bu (uz a ,uz b ) 



2 D 

P 



duu- 



yf*V$ ^U 2 ~k 2 -X 2 

G® (uz b , uz)G® (uz b , uz') cosh Uu 2 — k 2 — X 2 (t -if) . (2.40) 
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In this formula 

a J 

, K^\uzb) I v (uz)I v (uz') 
duu 



k*+\% ' ll b) (uz b ) ^ n 2 _ k 2 _ A 2 



x cosh 



u 2 -k 2 - \%t-t') (2.41) 



is the boundary part induced in the region z < z b by a single brane at z = z b when the brane 
z = z a is absent. Note that in the formulae given above the integration over angular part can 
be done by using the formula 



ikx / 



(2tt) 2 -A) (fc|x|)^- 

for a given function F(k). Combining two forms, formulae (|2.34j) and (|2.4Uj) . we see that the 
expressions for the Wightman function in the region z a < z < z\, is symmetric under the 
interchange a T^- b and l v <^ K v . Note that the expression for the Wightman function is not 
symmetric with respect to the interchange of the brane indices. The reason for this is that 
the boundaries have nonzero extrinsic curvature tensors and two sides of the boundaries are not 
equivalent. In particular, for the geometry of a single brane the VEVs are different for the regions 
on the left and on the right of the brane. For the case S = R Di one has ij)p(X) = (27r)~' D2 / 2 e jKX , 
where j3 = K is the ^-dimensional wave vector and replacing in the formulae above Yip ~^ 
(2tt)~ D2 J dK. we obtain the results for the AdSo+i bulk investigated in Ref. [SHI- Examining 
(|2.41|) tells one that the brane induced part in the Wightman function vanishes if one or both its 
arguments lie on the AdS boundary corresponding to z = 0. Note that this does not mean that 
the presence of the brane in the bulk has no consequences in the conformal field theory on the 
boundary, translated from the bulk physics by using the AdS/CFT correspondence. Here one 
should bear in mind that in the AdS/CFT correspondence the asymptotic behavior of the bulk 
modes <p a (x^, X 1 , z) is related to the corresponding boundary fields (po a (x^, X" 1 ) by the relation 
tpo a (x^, X 1 ) ~ lim z ^o z~ D / 2 ~ u ip a (x IJ ' , X 1 , z). It follows from here that the presence of the brane 
in the bulk leads to the additional contribution in the boundary two point function which is 
obtained from (|2.41j) as the limit Yim. z z t-+o(zz'}~ D ' 2 ~ v ((p(x)ip(x')) W. This limit is nonzero and 
gives nontrivial consequences in the boundary conformal field theory. It would be interesting 
to discuss the holographic interpretation of the brane degrees of freedom in terms of conformal 
field theory. However, such a discussion is beyond the scope of the present paper. 

Now let us consider the application of the results given above to the higher dimensional 
generalization of the Randall-Sundrum braneworld based on the bulk AdSo 1 +i x S. In the 
corresponding model y coordinate is compactified on an orbifold, S l jZ% of length I, with — I < 
y < I. The orbifold fixed points at y = and y = I are the locations of two D-dimensional 
branes. The corresponding line element has the form l|2.2j) with u{y) = p(y) = kr>\y\. The 
absolute value sign here leads to <5-type contributions to the corresponding Ricci scalar, located 
on the branes. As a result (5-terms appear in the equation for the eigenfunction f n (y)- Additional 
5-terms come from the surface action of the scalar field with mass parameters c\ and ci on the 
branes at y = and y = I respectively. The (5-terms in the equation for the function f n (y) in 
combination with the Z<i symmetry lead to the boundary conditions on this function. By the 
way similar to that for the usual Randall-Sundrum model (see, for instance, [5} I30| 139] ) it can 
be seen that for untwisted scalar the boundary conditions are Robin type with the coefficients 

^L = -h Cl +4D(k D ), ^ = -h-c 2 +4DCk D ). (2.43) 
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For twisted scalar field Dirichlet boundary conditions are obtained. Note that in the orbifolded 
version the integration in the normalization integral for the function f n (y) goes over two copies 
of the bulk manifold. This leads to the additional coefficient 1/2 in the expression for the nor- 
malization coefficient C n in (|2.24|) . Hence, the Wightman function in the higher dimensional 
Randall-Sundrum braneworld is given by formula 1)2.34(1 or equivalently by (|2.4U|) with an ad- 
ditional factor 1/2 and with Robin coefficients given by Eq. (|2,43|) , The one-loop effective 
potential and the problem of moduli stabilization in this model with zero mass parameters c,- 
are discussed in Ref. jlj. In particular, a scenario is proposed where supersymmetry is broken 
near the fundamental Planck scale, and the hierarchy between the electroweak and effective 
Planck scales is generated by a combination of redshift and large volume effects. The masses of 
the corresponding KK modes along £ are of the order of TeV. The phenomenology of this type 
braneworld models is discussed in Refs. [Hj ITU l IT5 | HI]. 

3 Vacuum polarization by a single brane 

In this section we will consider the VEV of the field square induced by a single brane located at 
y = a. As it has been shown in the previous section the Wightman function for this geometry 
is presented in the form 

<o|^(*Mz')|o) = M*Mz')> (0) + <^M*')> (a) , (3.i) 

where the boundary- free part on the left is defined by formula (|2.35j) . The brane induced 
part ((p(x)(p(x'))( a ^ is given by formula (|2.36l) in the region z > z a and by formula (|2.41|) with 
replacement Zb — > z a in the region z < z a . The VEV of the field square is obtained from the 
Wightman function taking the coincidence limit of the arguments. Similar to Eq. (|3.1|) this 
VEV is presented as the sum of boundary-free and boundary induced parts: 

<0|^|0) = ^ 2 > (0) + <^> (a) - (3-2) 

For the part without boundaries the direct substitution of the coincident arguments in the 
expression of the Wightman function leads to the divergent expression. For the regularization 
of this divergence the combination of the zeta function method for the sum over (3 and the 
dimensional regularization method for the integral over k can be used. By this way, taking the 
coincidence limit in the expression (|2.35j) and evaluating the k-integral, for the VEV of the field 
square in the boundary-free case one finds 

<^>(°) = ^|^r (i - ^) jHd„„c (i - ^,*;«) (3-3) 

where we have introduced the local spectral zeta function 

C(s,X;m) = y (3.4) 
p ( A | + m ) 

for the operator A^+^R^—m 2 . This corresponds to a scalar field with the mass m propagating 
on background of the manifold S. For £ = R D2 from (|3.3|) we obtain the standard result for the 
AdSo+i bulk. In this case the VEV of the field square does not depend on the coordinate z. 
For the bulk AdS x S, even in the case of homogeneous compact subspace S, this VEV depends 
on z. We explicitly illustrate this below in section [3 for the simple case £ = S . 

For points away from the brane the boundary induced part in the Wightman function is finite 
in the coincidence limit and we can directly put x = x'. Introducing a new integration variable 
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v = yjv? — k 2 — A^,, transforming to the polar coordinates in the plane (v,k) and integrating 
over angular part, the following formula can be derived 



oo poo 

dk 



k^+xi ^J u 2 -k 2 -\} 2r (Si) 

By using this formula and Eq. (|2.3fij) . the boundary induced VEV for the field square in the 
region z > z a is presented in the form 

<^> (a) =ElM*)lV>?, (3-6) 

P 

where the contribution of a given KK mode along X is determined by formula 

I N k D ~ 1 7 D D, T^(ll7 ) 

P 2 D i- 1 Tr^r(^)Jx Ki a) (uz a ) 

The corresponding formula in the region z < z a is obtained from Eq. H2.41|) by a similar way 
and differs from Eq. ()3.7|) by the replacements I u <^ K u . Using the properties of the Bessel 
modified functions, it can be seen that (ip 2 )^ < for Dirichlet boundary condition (B a = 0) 
and {(p 2 )^ > for the ratio of the coefficients in the interval — v < A a /B a < v in both regions 
z < z a and z > z a . In the general case, the VEV may change the sign as a function of z. In the 
models with a homogeneous internal space E, the sum over (3 in Eq. (|3.6|) and, hence, the VEV 
of the field square do not depend on the coordinates X 1 in this space. 

For the comparison with the case of the bulk spacetime AdSo 1 +i when the internal space 
is absent, it is useful in addition to the VEV (|3.6j) to consider the VEV integrated over the 
subspace X: 

<Al gra ted = / d^X^^e-^ = e -^*Xy>M (3.8) 

Comparing this integrated VEV with the corresponding formula from Ref. jHS], we see that the 
contribution of the zero KK mode (Xp = 0) in Eq. (|3.8j) differs from the VEV of the field square 
in the bulk AdSox+i by the order of the modified Bessel functions: for the latter case v — * v\ 
with v\ defined by Eq. (|2.19j) with the replacement D — ► D±. Note that for £ < C-D+-D1+1 one 
has v > v\. In particular, this is the case for minimally and conformally coupled scalar fields. 
If the internal space X is a one-parameter manifold of size L then one has ~ 1/L. By taking 
into account that from the normalization condition it follows that the function ij)a{X") contains 
the factor L~~ D2 ^ 2 , we see that ((p 2 )^ is a function on z/z a and L/z a . The first ratio is related 
to the proper distance from the brane by the equation 

z / Za = e kD{y-a)_ ( 39) 

Hence, for a given size L of the internal space, the VEV of the field square in addition to 
the proper distance from the brane, also depends on the absolute position of the brane in the 
bulk. Note that in the case of AdS bulk the corresponding quantity depends only on the proper 
distance from the brane [39] . To discuss the physics from the point of view of an observer 
residing on the brane y = a, it is convenient to introduce rescaled coordinates x ,fl = e~ k ° a x^. 
For this observer the physical size of the subspace S is L a = Le~ k ° a and the corresponding KK 
masses are rescaled by the warp factor, = \pe k ° a . Now we see that the VEV induced by 
a single brane is a function of the proper distance from the brane and on the ratio L a /(l//cc>) 
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of the physical size for the internal space E for an observer residing on the brane to the AdS 
curvature radius. 

As a partial check for derived formulae let us consider the limit kjj — * 0. This corresponds 
to a boundary on the bulk R^ 1 ' 1 ' x E. For kr> — ► 0, from Q2.19|) we see that the order u of the 
cylindrical functions is large. Introducing in ()3.7j) the new integration variable v = u/u, we can 
replace the Bessel modified functions by their uniform asymptotic expansions for large values of 
the order. To the leading order this gives: 



P 

Di P - 2u \y-<A 
_du(u 2 -m 2 - Xl)^- 1 -— — . (3.10) 



v /m2 +A 2 P C a {u) 



Here and below we use the notations 



d j( u ) = lT F7w^' 3 = a,b, (3.11) 



Aj + nO)B 



u 



and is defined after formula (|2.38j) . 

For the points on the brane the VEV given by formula H3.7|) diverges due to the divergence 
of the u-integral in the upper limit. The surface divergences in the renormalized VEVs of the 
local physical observables result from the idealization of the boundaries as perfectly smooth 
surfaces which are perfect reflectors at all frequencies, and are well known in quantum field 
theory with boundaries (see, for instance, [5H1 031 ES] ) • ^ seems plausible that such effects as 
surface roughness, or the microstructure of the boundary on small scales can introduce a physical 
cutoff needed to produce finite values of surface quantities. In brane models the imperfectness 
would come from the quantum gravity effects on the Planck scale. An alternative mechanism for 
introducing a cutoff which removes singular behavior on boundaries is to allow the position of the 
boundary to undergo quantum fluctuations [HUJ. Such fluctuations smear out the contribution 
of the high frequency modes without the need to introduce an explicit high frequency cutoff. 
In order to find the renormalized value of the field square on the brane many regularization 
techniques are available nowadays. In particular, the generalized zeta function method is in 
general very powerful to give physical meaning to the divergent quantities. This method has 
been used in Ref. [IS] for the evaluation of the field square on the brane in the model with 
AdSn+i bulk. Another possibility is to remove from (|3.7|) the terms diverging in the limit 
z —* z a . The relation between various procedures for the evaluation of the renormalized VEV of 
the field square on the brane in the thin wall approximation is recently discussed in Ref. |49| . 
In particular, on the specific example it has been shown that different methods give the same 
result up to finite renormalization of the surface mass and extrinsic curvature counter-terms. 

From the divergent behavior of the vacuum fluctuations on the brane it follows that for 
the points near the brane the main contribution into the integral comes from large values of it. 
Assuming \p\z — z a \ <C 1, holy — a\ <C 1 and replacing the Bessel modified functions by their 
asymptotic expansions for large values of the argument, to the leading order for the contribution 
of the mode with a given (3 one finds 

,<■.,,_ (W>-«(g,) r /A_n 

p {^)^\z- Za \ D ^ V 2 ; 



where the notation 



K{B a ) = 25 Ba0 ~ 1 (3.13) 
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is introduced. Note that in this limit the proper distance from the brane is much smaller 
compared with the AdS curvature radius and the proper size L a of the internal manifold E. 

Now let us consider the contribution of the large KK masses along the internal subspace 
S, assuming that z a Xp,zX/3 ^> 1. This corresponds to the KK masses \$ much larger than 
the AdS energy scale hp. In particular, for L a <C 1/&d these conditions are satisfied for all 
modes along £ with nonzero KK masses. Introducing a new integration variable v = u/Xp, we 
can replace the Bessel modified functions by their asymptotic expressions for large values of the 
argument. For the contribution of the mode with a given (3 one finds 

(V% a) *- ( T — / du^-X^^-y. (3.14) 
' (47r)^r (^l) JXp c a {uz a ) 

where we have introduced the notation 

. . A, - n^BjU 
c i W = — m"^' 3 = a ' 6 > ( 3 - 15 ) 



and is defined after formula 1)2.38)1 . Additionally assuming that z a Xn 3> A a /B a or £> a = 0, 
one has 



2 Dl 7T 2 



Xl 



K £l -i(2Xp\z-z a \). (3.16) 



From this formula it follows that for Xp\z — z a \ 3> 1 (the observation point is not too close to 
the brane) the contribution of KK modes with large is exponentially suppressed. 

For small values z a , z a <C z,l/Xp, using the asymptotic formulae for the Bessel modified 
functions for small values of the argument, we obtain the formula 

tf)f - ~ \ *° Za / duu**V - Xl)^Kl{uz). (3.17) 



This limit corresponds to large distances from the brane, kz>(y — a) 3> 1 and small KK masses 
\p compared with the AdS energy scale, A^ a 



Xg compared with the AdS energy scale, aL <C &d- In particular, from Eq. (|3.17|) it follows 



that for fixed values ko, y, A/3 the brane induced VEV in the region z > z a vanishes as z 2 " when 
the brane position tends to the AdS boundary, z a — > 0. Formula 1)3.1 7|) is further simplified for 
two subcases. In the case z a <C A^ 1 <C z or equivalently k£>e~ kD<yy ~ a ^ <C Xy <C ko, the main 
contribution into the integral comes from the lower limit ant to the leading order we find 

£1 

{ip2) {a) ra 1 y^p) 2 e -2zx^ (3 _ 18) 



k D D ~\^{ZqX^{zX^-\ _ 2zXh 



with the exponential suppression of the boundary induced VEV. In the limit z a <C z -C A^ 1 

(small KK masses, x f 1 -C k£ie^ kD ^ y ^ a ^) to the leading order we can put in the lower limit of 
the integral. Then the integral is evaluated by the standard formula for the integrals involving 
the square of the MacDonald function and one finds (see, for instance, |61j ) 

(^- A'-f r{ l + 11»Sf + : ] (zf ■ (3.1.) 



In particular, this formula is valid for the zero mode. 
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In the limit z <C z a , again using the asymptotics of the modified Bessel functions for 

small values of the argument, we obtain the formula 



\uz, 



a J 



This limit corresponds to large proper distances from the brane k£>{a — y) S> 1 in the region 
y < a. As we see for fixed values ko, z a , \r the brane induced VEV vanishes as z D+2v when the 
observation point tends to the AdS boundary. As it has been explained in the previous section 
[see paragraph after formula (|2.42|) ]. this does not mean that the presence of the brane in the 
bulk does not contribute to the corresponding VEV in the boundary conformal field theory 
within the framework of AdS/CFT duality. To further simplify the formula (j3.20j) we consider 
two subcases. For the case z <C A^ 1 -C z a corresponding to ^-oiy-a) ^ k D jy^^> <^ l ; to the 
leading order we obtain the following formula 



fcg- 1 z^c a (z Q A /3 )e- 2 ^ A ^ M ^ 

2£i+2*-i 7r ^- 1 r 2 0/ + i) 



Ex 
z \ - 

Z a, 



In the opposite limit for z a , z <C z a <C Ag 1 , which corresponds to small KK masses, Li- 
the lower limit in the integral can be replaced by and to the leading order the contribution of 
the mode with a given f3 does not depend on A^: 

^r^r^ + i) W io ^(n) 

Next we consider the limit z\p S> 1 when ZaA^ is fixed, corresponding to e fcl3( ^~ a * ) 3> kjj/X^, 
y > a. The main contribution to the integral in Eq. (|3.7|) comes from the lower limit and to the 
leading order one finds 

^ „ - (foz)g ::f:;^ e-^, (3 , 3) 

with the exponentially suppressed VEV. In particular, the contribution of the nonzero KK modes 
along S exponentially vanishes when the observation point tends to the AdS horizon, z — » cxd. 
Note that for the zero mode (A/3 = 0) the brane induced VEV near the AdS horizon behaves as 
z D 2 -2u ^ gee f ormu i a (j3.i9jQ . In the purely AdS bulk (D2 = 0) this VEV vanishes on the horizon 
for v > 0. For an internal spaces with D2 > 2u the VEV diverges on the horizon. Note that 
for a conformally coupled massless scalar and D2 = 1 the boundary induced VEV takes nonzero 
finite value on the horizon. 

In the limit when the brane position tends to the AdS horizon, z a — > 00, for massive KK 
modes along £ the main contribution into the VEV of the field square in the region z < z a 
comes from the lower limit of the u-integral. To the leading order we find 



W « / g Z\ P \ > (3-24) 

2°1 IT^Za 2 



and the VEV is exponentially small. As it can be seen from general formula, for the zero mode 
in the same limit the VEV vanishes as z~ Dl ~ 2u . 
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It is also of interest to consider the behavior of the VEV for large values of the AdS energy 
scale ko corresponding to strong gravitational fields. When the values of the other parameters 
and the coordinates a and y are fixed, for nonzero KK modes one has X^z a , \pz 3> 1. Addi- 
tionally assuming \p\z — z a \ S> 1, from Eq. (|3.16|) we see that (<p 2 )a behaves as (koz) D ~ l \z — 
Za\~ Dl ^ 2 ^xp[—2\p\z — z a \] and is exponentially small. For the zero KK mode from the general 
formulae we can see that the corresponding contribution to the VEV of the field square behaves 
as k^~ x e D ^ k ^ exp[2vk D (a-y)] in the region y> a and like k^~ x e D ^y exp[(D 1 +2u)k D (y-a)} 
in the region y < a. Note that the corresponding factors in the VEVs integrated over the in- 
ternal subspace (see Eq. (|3.8jl ) contain an additional factor e~ D2k ° y coming from the volume 
element and are exponentially suppressed everywhere in the parameter space. Hence, the strong 
gravitational field suppresses the vacuum fluctuations. The similar behavior for the boundary 
induced quantum effects in the gravitational field of the global monopole is described in Ref. |62j . 



4 VEV in two branes geometry 

The branes divide the space into three distinct sections: < z < z a , z a < z < z&, and z > Z\>. 
In general, there is no reason for the inverse curvature radius to be the same in these three 
sections, as branes may separate different phases of theory. The VEVs in the first and last regions 
are the same as in the corresponding geometry of a single brane described in the previous section. 
In this section we consider the VEV for the field square in the region between two branes. Note 
that in the braneworld scenario with two branes based on the orbifolded version of the model 
this region is employed only. Taking the coincidence limit in the corresponding formulae for the 
Wightman function and using the integration formula (|3.5j) j for the VEV of the field square we 
obtain two equivalent representations 

«V|o> = <^)(°) + <^>w k -K^— T £hfc(*)l 2 

20x-%^r(SL) V 



roc D 

x / duu(u 2 - \ 2 p)^~ 1 Q, jv (uz a ,uz b )G^ 2 (uz j ,uz), (4.1) 
J\ g 

with j = a, b. The last term on the right of this formula is finite on the plate at z = Zj and 
diverges for the points on the brane z = Zj>, f ^ j. These divergences are the same as those for 
a single brane at z = zy . As a result if we write the VEV of the field square in the form 

j=a,b 

then the interference part (y> 2 )( a6 ) is finite on both plates. In the case of one parameter manifold 
£ with the size L and for a given hp, the VEV in Eq. ()4.1j) is a function on Zb/z a , L/z a , and 
z/z a . Note that for the points on the brane z = Zj in the last term on the right (|4,lj) one has 
Gu \uzj,uzj) = —Bj and this term vanishes for Dirichlet boundary conditions. 

By the same way, as in the single brane case, it can be seen that in the limit ko — ¥ the 
result for the geometry of two parallel Robin plates in the bulk R> Dl '^ x S is obtained: 

2\(o6) _ ( 47r )~ 



(u 2 - m 2 - A 2 ,)*-! 

du — , . n „ ; 

1 



y/^f c a {u)c h (u)e 2 < b - a ) 



e — 2«[s/— jI 



Cn (it) 

j=a,b JV ; 



(4.3) 
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with the notation Cj(u) from Eq. (|3.11l) . 

For large KK masses along S, z a Xp S> 1, for the contribution of the mode with a given /3 
one finds 



(k D z 



da 



A 



(4vr)^ I r (Sl) c a {uz a )c b {uz b )e 2 < z b z °) - 1 



£ 

j=a,b 



,—2u\z—Zj | 
Cj(uZj) 



(4.4) 



where (f 2 )^^ is defined by the relation similar to ()3.6|) . Additionally assuming the conditions 
A / 3(zfo — z a ) 3> 1 and Xp\z — Zj\ 3> 1, for this quantity we obtain the formula 



2\(a6) _ (kpz) 



D-l 



-2\p{z b -z a ) 



V 



(4.5) 



In particular, for L a <C 1/£;.d this formula takes place for all KK modes along £ with nonzero 
masses. 

Now let us consider limiting cases when the general formula for the interference part (ip 2 )( ab ) 
can be simplified. First of all assume that \pz b S> 1 and z a Xp < 1. This limit corresponds to 
large interbrane distances compared with the AdS curvature radius k~j^ and is realized in the 
braneworld scenarios for the solution of the hierarchy problem. In this limit the main contribu- 
tion into the u-integral comes from the region near the lower limit and for the interference part 
to the leading order we have 



k^- 1 z D X 2 e~ 2X ^ 



2 D ^^- 1 z b 2 

li a) (Zg\p) 

Ki a \z a Xp) 



c b (\pz b ) 



■K v {\pz) 



Z a X(l) 



Kl a \z a Xp] 



K v {X p z) - 2I v (\ p z) 



(4.6) 



for the nonzero KK modes along S. In the limit z a , z <C X^ 1 <C z b the corresponding formula 
takes the form 



k^z^jzaX^jzX^e- 2 ^ 
2 D ^-^- 1 c a (u)c b (z b X p )T 2 (u + 1) U 



2c a (v) 



2ul 



(4.7) 



For Xpz 3> 1 and z a Xa < 1, by using asymptotic formulae for the Bessel modified function one 
finds 



; 2\(ab 



k^z^jzX^- 
2 D ^~ l c b {z b X p ) 



li a \zgXp) 

Ki a \z a X p ) 



e -2\ fj (z b -z) 
Xf3Z b C b (X[3Z b ) 



(4.8) 



And finally, in the limit z a <C z b <C A^ 1 to the leading order we can put instead of Xp in the 
lower limit of the integral over u and by using the asymptotic formulae for the Bessel modified 
functions for small values of the argument, it can be seen that {<-p 2 )^ ~ {z a / z b ) 2u g(z/ z b ). If in 
addition one has z <C z b the following formula is obtained 

'Z^izjz^iz/Z,)^ 



I 2\{ab 



k D D l ^l 



2v 



2Cg{v) 



du u 



D 



1+21/-1 \u) 

ll b \u) 



(4.9) 
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with the exponential suppression of the interference part. As we see for these values of the 
parameters the interference part in the VEV of the field square is mainly located near the brane 
at z = Z\). In particular, from formula (|4,9jl it follows that for the zero KK mode along X 
the interference part in the VEV of the field square vanishes as z^ Dl ~ 2u when the right brane 
tends to the AdS horizon. Note that in the same limit the contribution of a given massive 
KK mode vanishes as e~ 2X/3Zb / zj^ 1 ^ 2 (see formula (|4.6|) ). When the left brane tends to the AdS 
boundary, z a — ► 0, the interference part vanishes as z 2u . The behavior of the interference part 
in the VEV of the field square for large values of the AdS energy scale kn directly follows from 
formula (|4.5jl in the case of nonzero KK masses and from formula (|4.9j) in the case of the zero 
mode. Note that the corresponding contributions to the VEV of the field square integrated over 
the internal space contain additional factor e~ D2k ° y coming from the volume element. In the 
scenario considered in Ref. L < 1/ku, z a ~ L, and z^jz a 3> !• For these values of the 
parameters the contribution into the interference part for the VEV of field square to the leading 
order is given by Eq. ()4.6|) for nonzero KK modes and {^p 2 }^^ ~ (z a /zi ) ) 2l/ for the zero mode. 
If in addition the observation point is far from the brane at z = z^, the zero mode contribution 
is determined by formula (|4,9[) with an additional suppression factor (z/zb) Dl . 



5 An example 

In the discussion above we have considered the general case of the internal space. Here we take 
a simple example with X = S . In this case the bulk corresponds to the AdSn+i spacetime 
with one compactified dimension X. The length of this dimension we will denote by L. The 
corresponding normalized eigenfunctions are as follows 

MX) = -±=e 2 ^ x ' L , \p = ^-\(3\, P = 0,±l,±2,.... (5.1) 

First of all we consider the Wightman function for the bulk without boundaries given by formula 
(|2.35|) . We apply to the sum over (3 the Abel-Plana summation formula 



£'/(/*) 



f3=0 



oo 



■*)+'• i ^ ^t/ZV^ (5 - 2) 







where the prime on the sum sign means that the summand (3 = should be taken with the 
weight 1/2. Now it can be easily seen that the term in the Wightman function with the first 
integral on the right of formula (|5.2|) corresponds to the Wightman function ((p(x)ip(x'))^ > dSD+i 
for the scalar field in the bulk AdSp+i- The latter is well investigated in literature. As a result 
for the Wightman function in the bulk AdSr> x S 1 one finds 

M*M*')> = M*M*))i°k +1 + k \n-^TS I dke*^ 
dv vJ v {vz)J v (vz) 







OCj 







, cosh[u(t - t')] cosh[(X - X'Wu 2 + v 2 + k 2 } . . 

du , ; „ = = . (5.3 

Vu 2 + v 2 + k 2 eWt^+t^+fc* _ i v ; 



To evaluate the corresponding VEV for the field square we take in this formula the coincidence 
limit. In this limit all divergences are contained in the part { L p(x)ip(x'))^j lSD+i , and the additional 
part on the right of formula ()5.3|) induced by the compactness of a single dimension leads to 
the finite result. To evaluate the corresponding integrals in the coincidence limit, first we 
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introduce polar coordinates (r, 9) in the plane (k,u). After the evaluation the simple ^-integral, 
we introduce the polar coordinates (v, <j>) in the plane (x,r). By using the formula 61 

1 , 1 o, , , , /«\ & 1F2 (v+ 1/2: f + 1 + (3,2v + l;-v 2 ) , . 



the integral is expressed through the hypergeometric function 1F2, and one finds 

2,(0) , 2k D D ~\z/2L) D +^ 



\Y I \Y IAdSo+i^ ±!=±tV , i\r/ , D+V 



vr— I> + l)r {u+S±k) 

roc u D+2u-\ / 1 D + 1 Z 2 U 2\ 

X J Q du e u_ 1 lF2 ^ + 2 ;iy + ^~' 2l/ + 1; ~T^j- (5 - 5) 

As we see, unlike to the case of AdSn+i bulk, the VEV for the field square in the bulk AdSo x S 1 
is a function on z. Note that the second term on the right of formula (|5.5|) is always positive. For 
small values of the ratio z/L, to the leading order the integral is equal to F(D + 2v)(r(D + 2u), 
with Cr( x ) being the Riemann zeta function, and this term behaves as (z / L) D+2u . For large 
values of z/L it can be seen that to the leading order one has 

where we have taken into account that the quantity (^P/AdSni does n °t depend on z. Hence, 
near the AdS horizon the VEV of the field square is dominated by the second term on the 
right of Eq. (|5.5|) . Note that in the limit — > the VEV ( [ P 2 )AdS D+1 varns hes, and taking 
into account that zko — ► 1, from (|5.fij) we obtain the standard result for the VEV in the bulk 

R (D-1,1) x 51. 

Brane induced VEVs of the field square for the geometry under consideration are obtained 
from general formulae (|3.7|) . (|4.1j) by the replacements 

52\M*)\ 2 - '» D^D- 1 - (5-7) 

13 f3=0 

As it has been mentioned in section |21 for D2 = 1 and conformally coupled massless scalar 
field {y = 1/2) the single brane induced VEV is finite on the AdS horizon. The corresponding 
limiting value is obtained from (J3.19|) : 

, av(.). , _ k^T(D/2)z a 
2 D -H^c a (l/2)L 

For the general case, the VEV in this limit behaves as z 1-2 ^. In figure ^ we have plotted the 
dependence of the single brane induced VEV (<p 2 )( a ' on the ratio z/z a for D = 5 minimally (left 
panel) and conformally (right panel) coupled massless scalar fields with A a /B a = —0.4. The 
full and dashed curves correspond to the values L/z a = 1 and L/z a = 0.2 respectively. As it has 
been shown in section 03 in the limit z — > the VEV behaves as z D+2u . For the points near the 
AdS horizon, z — ► 00, the main contribution comes from the zero mode and the VEV behaves 
as z 1-2 ^. On the brane surface the VEV diverges like \z — z a \ 2 ~ D . 

In figure EJ we present the graphs for the interference part in the VEV of the field square 
as a function on z/z a for D = 5 minimally (left panel) and conformally (right panel) coupled 
massless scalars with Aj/Bj = —0.4, j = a, b. The graphs are plotted for the interbrane distance 
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Figure 1: Single brane induced VEV of the field square as a function on zj z a for D = 5 minimally 
(left panel) and conformally (right panel) coupled massless scalars with A a /B a = —0.4. The full 
and dashed curves correspond to the values L/z a = 1 and L/z a = 0.2 respectively. 




Figure 2: Interference part in the VEV of the field square as a function on z/z a for D = 5 
minimally (left panel) and conformally (right panel) coupled massless scalars with Aj/Bj = —0.4, 
j = a,b. The full and dashed curves correspond to the values L/z a = 1 and L/z a = 0.2 
respectively. The interbrane distance corresponds to the value Zb/z a = 3. 
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corresponding to z^/ z a = 3. Note that in accordance with the general discussion in the previous 
section one has (v9 2 )i=l a < ((p 2 ){=l b . 

It is not difficult to generalize the corresponding results for the case of the internal space 
E = S° 2 with radius Rq and curvature scalar Rr^ = — D 2 {D2 — Now the eigenfunctions 

tpp(X) are expressed in terms of spherical harmonics of degree I, I = 0, 1, 2, The VEVs for 

the field square are obtained from the general formulae in sections |3] and 21 by the replacements 



Z7T ^ JTq i=0 

A/3 - 



+ D 2 -l) +CD 2 (D 2 -l), (5.10) 

"0 



where the factor under the sum sign on the right in Eq. (|5.9|) is the degeneracy of each angular 
mode with a given I. 



6 Conclusion 

From the point of view of embedding the braneworld model into a more fundamental theory, such 
as string/M-theory, one may expect that a more complete version of this scenario must admit 
the presence of additional extra dimensions compactified on a manifold E. In the present pa- 
per we have extended the previous work describing the local vacuum effects in the braneworlds 
with the AdS bulk on a higher dimensional brane models which combine both the compact 
and warped geometries. This problem is also of separate interest as an example with gravita- 
tional, topological, and boundary polarizations of the vacuum, where one-loop calculations can 
be performed in closed form. We have investigated the Wightman function and the vacuum 
expectation value of the field square for a scalar field with an arbitrary curvature coupling pa- 
rameter satisfying Robin boundary conditions on two parallel branes in AcLSdi+i x S spacetime. 
The KK modes corresponding to the radial direction are zeros of a combination of the cylinder 
functions. The application of the generalized Abel-Plana formula to the corresponding mode 
sum in the expression of the Wightman function allowed us to extract the boundary-free part 
and to present the brane induced parts in terms of integrals exponentially convergent in the 
coincidence limit of the arguments. In the region between two branes the Wightman function is 
presented in two equivalent forms given by Eqs. 1)2.341) and (|2.4U[) • The first term on the right 
of these formulae is the Wightman function for the bulk without boundaries. The second one is 
induced by a single brane and the third term is due to the presence of the second brane. Further 
we give an application of our results for two branes case to the higher dimensional version of 
the Randall- Sundrum braneworld with arbitrary mass terms on the branes. For the untwisted 
scalar the Robin coefficients are expressed through these mass terms and the curvature coupling 
parameter by formulae 1)2.43)) . For the twisted scalar Dirichlet boundary conditions are obtained 
on both branes. 

The expectation values for the field square are obtained from the Wightman function taking 
the coincidence limit. For the case of a single brane geometry this leads to formula ()3.7I) for the 
region z > z a . The corresponding formula in the region z < z a is obtained from Eq. 1)3. 7jl by 
replacements I u T± K u . For a one parameter manifold E, the VEV of the field square induced 
by a single brane is a function on the proper distance of the observation point from the brane 
and on the ratio of the physical size of E (from the viewpoint of an observer residing on the 
brane) to the AdS curvature radius. As a partial check of our formulae, we have shown that in 
the limit when the AdS radius goes to infinity the result for the brane in the bulk R^ 1 ' 1 ) x E is 
obtained. On the boundary the VEV for the field square diverges. In the limit when the proper 
distance from the brane is much smaller compared with the AdS curvature radius and proper 
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size L a of the internal manifold, the leading term of the corresponding asymptotic expansion 
for a contribution of a given KK mode along £ is given by expression (|3.12|) . This term does 
not depend on the Robin coefficient and has different signs for Dirichlet and non-Dirichlet 
scalars. The coefficients for the subleading asymptotic terms will depend on the bulk curvature 
radius, Robin coefficient and on the mass of the field. In the limit when the KK masses aL 
are much larger compared to the AdS energy scale kjy, the contribution of a given mode with 
fixed (3 to the VEV of the field square is given by formula ()3.16|) . From this formula it follows 
that if additionally one has A^|y — a\ 3> 1, the contribution of the corresponding KK modes is 
exponentially suppressed. In the region y > a, for large distances from the brane, ko{y — a) S> 1, 
and small KK masses A^ compared with AdS energy scale, Al <C k^, the contribution of a 
given KK mode along the internal space £ is approximated by formula (j3,17|) . In both subcases 
kne~ kD ^ y ~ a ^ <C and ki)e~ kD ^ y ~ a ^ 3> Ao we have exponentially suppressed VEVs by the 

factors exp [-2e kD{ -y-^\ ( p ] /k D ] and exp[2v(a — y)] in the first and second cases respectively. 
For large distances from the brane in the region y < a, the contribution of a given KK mode 
is presented in the form (|3.2UI) . This formula is further simplified in the subcases of large 
and small KK masses A^. In the case of large mass, ^-oiy-a) <g fcp/A'g <C 1, one has an 

exponential suppression by the factor exp[— 2A^//cd]. In the limit of small masses, \f « k D, 
the suppression is by the factor exp[ku(y— a)(D\+2v)\. On the AdS horizon (z — > oo), the brane 
induced VEV is exponentially small for nonzero KK masses along £ and behaves as z D ' 2 ~ 2u for 
the zero mode. On the AdS boundary, corresponding to z = 0, the brane induced VEV vanishes 
as z D+2u . When the brane position tends to the AdS boundary, z a — > 0, the brane induced VEV 
in the region z > z a vanishes as z 2u . In the limit when the brane tends to the AdS horizon, 

D /2 

z a — > oo, the VEV in the region z < z a behaves as exp[— 2\pz a ]/z a for massive KK modes 
along E, and vanishes as z~ Dl ~ 2u for the zero mass mode. For large values of the AdS energy 
scale ko corresponding to strong gravitational fields, the VEVs integrated over the internal 
space (see formula (j3.8[l ) are exponentially suppressed everywhere in the parameter space. 

For two branes geometry, by using the corresponding expression for the Wightman function, 
VEV for the field square in the region between the branes is presented as a sum of boundary- 
free, single brane induced and interference parts, Eq. (|4.2j) . The latter is regular everywhere 
including the points on the branes. In the limit kp — > the result for the geometry of two 
branes in the bulk R^ Dl ^ x S is obtained. For large KK masses along S, 3> kp, the 
contribution of the KK mode with a given (5 is determined by formula (|4,4j) . This formula is 
further simplified when the branes are not too close to each other and the observation point is 
sufficiently far from the branes, \p{zt> — z a ) 1 and \p(z — zj) 3> 1, j = a,b. In this case 
we obtain formula ()4.5|) with exponentially small interference part. For the modes with masses 
in the range z«A^ <C zj,, the contribution to the VEV of the field square is determined by 
formula ()4.7j) with two small factors, e _2A ' 52b and (z / ' z^) Dx l 2 . For small values of KK masses, 
Za *C zi, A7 , and large interbrane distances, to the leading order the contribution of a 
given KK mode does not depend on the KK mass and the corresponding VEV is exponentially 
suppressed by the factor exp[2vkD(a — b)]. If in addition the observation point is far from the 
brane at y = b, additional suppression factor exp[Z)i&;£)(y — b)] appears and the interference part 
is mainly located near the brane y = b. When the right brane position tends to the AdS horizon, 
Zb — > oo, the interference part in the VEV of the field square vanishes as z^ Dl ^ 2u for the zero 

KK mode along S. For a given massive KK mode this quantity behaves as exp[— 2\pz ]/ z® 1 ^ 2 . 
When the left brane tends to the AdS boundary, the interference part vanishes as z 2u . 

As an illustration of general results, in section 03 we consider an example of the internal space 
S = S . Using the Abel-Plana formula, the Wightman function for the boundary-free geometry 
is presented in the form (|5..3ft . where the first term on the right is the corresponding function 
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for the AdSp+i bulk. The VEV of the field square is obtained in the coincidence limit and has 
the form ()5.5|) . Unlike to the case of AdSp+i bulk, this VEV for the bulk AdSo x S 1 depends 
on z coordinate and behaves as (z/L) D+2u for small values of the ratio z/L and as (z/L) ^ 1 
for large values of this ratio. In figures ^ an d 121 we have presented single brane induced and 
interference parts in the VEV of the field square for minimally and conformally coupled scalars 
in the bulk geometry AdS§ x S 1 as functions on the distance from the brane at y = a. We also 
discuss the generalization of the corresponding results for the internal space S D2 . TheVEVs m 
this case are obtained from the general formulae in sections |3] and |1] by the replacements (|5.9f) 
and (|5.10[) . In particular, from the point of view of embedding the Randall-Sundrum model into 
the string theory, the case of the bulk AdS§ x S" 5 is of special interest. This spacetime also plays 
an important role in the discussions of the holographic principle. 

Note that in this paper we have considered boundary induced vacuum densities which are 
finite away from the boundaries. We expect that similar results would be obtained in the model 
where instead of externally imposed boundary condition the fluctuating field is coupled to a 
smooth background potential that implements the boundary condition in a certain limit |63| . 
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